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Figure 1: exterior domain with artificial interface T 

1 Introduction 

Tlie main focus of our investigations is to suggest a method of deriving guaranteed and 
computable upper bounds of the difference between the exact solution u of an elliptic 
exterior domain boundary value problem and any approximation from the corresponding 
energy space. We discuss the method with the paradigm of the prototypical elliptic 
problem 

-divAVu = f inn, (1.1) 

= g on 7 := dQ. (1.2) 

We assume that Q C with > 1 is an exterior domain, i.e. \ is compact, with 
Lipschitz continuous boundary 7 (see Figure [T|. 

Throughout this paper we will use the weighted Lebesgue function spaces 

ll{n) := I pV e l\n)} , seR. 

Here p := (1 + r^)^/^ and r{x) := \x\ denotes the radius vector. L^(fi) is a Hilbert space 
equipped with the scalar product 

Jn 

where and ip belong to L^(f2) and A is Lebesgue's measure. We denote the corresponding 
norms by = \\p^(p\\q. If s = then L^(r2) coincides with the usual Lebesgue space 



A Posteriori Error Estimates in Exterior Domains 



3 



L^(n). For the sake of simplicity we keep the same notation for spaces of vector- valued 
functions. Moreover, we introduce the weighted Sobolev space 

which is a Hilbert space as well with respect to the scalar product 

(v^, ij) ^ {if, + {v^, Vi^)^ . 

o o 

By H^]^(fi) we denote the closure of C°°(fi), the space of compactly supported smooth 
test functions, in the norm of Hi]^(i7). Whenever we consider Sobolev spaces of bounded 
domains we use the usual unweighted L^-scalar products and -norms. 

For dimensions > 3 the solution theory for the problem (1.1)-(1.2) is based on 
the weighted Poincare/Friedrich estimate (see Corollary 16 (i) and Remark 17 of the 
appendix) 



M-in < 



;i.3) 



the Lax-Milgram theorem and, if needed, an adequate extension operator for the boundary 
data. Let be some function in H^^(fi) satisfying the boundary condition (1.2). The 



weak solution u G H|'Li(i7) + C H|'Li(i7) of (1.1)-(1.2) is then defined by the variational 
formulation 



(1.4) 



By (1.3) the left hand side of (1.4) is a strongly coercitive sesqui-linear form over H|'Li(r2) 
provided that the real-matrix- valued function A is measurable, bounded a.e., symmetric 
and uniformly strongly elliptic, i.e. 



3ca>0 y^eR^ VxG^] ca\^\' < A{x)^ ■ t 



;i.5) 



If / G L]^(f2) then by the Cauchy-Scharz inequality the right hand side of (1.4) is a linear 



and continuous functional over Hi]^(n). Thus, under these assumptions the problem (1.4) 

o 

is uniquely solvable in Hi]^(f2) -|- by Lax-Milgram's theorem. 

If = 1, 2 one can apply the same arguments with the difference that (1.3) has to be 
modified. For N = 1 and, for example, Q C M+ we have by Corollary 16 (iii) and Remark 



llv'ILi n < 2 iv^lln V^GHii(fi). 



;i.6) 



Hence, we get the same solution theory with tiny restrictions on Q, which easily can be 
removed by a translation. For N = 2 the singularities are stronger and additionally we 
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have to utilize logarithmic terms. By Corollary [16] (ii) and Remark 17 we have for domains 
Q C M^, such that the complement \ i7 contains the unit ball, 

||Mrlnr)||^<2||V^||^ \/cpeH\U^), (1.7) 

where 

:= I ^/(rlnr),V^GL2(n)} 
is a Hilbert space equipped with the natural scalar product 

(<^,^) ^ (<^/(rlnr), V/(rlnr))f^ + {V(p,V^)^ 

o o 

and again H^^^ in(^) denotes the closure of C°°(fi) in the norm of H^^^ in(^)- Consequently, 
we obtain for all / with rlnr/ G L^(r2) and all u-y in Hi.^jjj(fi) satisfying the boundary 



condition (1.2) a unique solution u belonging to H^^^ jj^(r2) + u^. 
We summarize the results in the following 

Theorem 1 Let N > 3 as well as f G L^(fi) and E H\^(n) satisfying the boundary 
condition (1.2). Then the exterior boundary value problem (1.1)-(1.2) is uniquely weakly 

o 

solvable in Hi.j^(i7) +u^. The solution operator is continuous. 

From the above discussion, it is clear that for = 1, 2 the existence of weak solutions 
in suitable spaces can also be proved. 

Remark 2 The boundary data g and its extension u-y can be described in more detail. In 
the bounded domain case it is well known that there exists a bounded linear trace operator 
and a corresponding bounded linear extension operator (right inverse) mapping to 
H-^/^(7) and vice verse. Hence, by restriction we get a bounded linear trace operator 

and by extension and applying an obvious cutting technique we obtain a bounded linear 
extension operator 

E : H^/2(^) ^ |_|^^(^) 

for our exterior domain Q, which even maps to functions with (arbitrarily thin) compact 
support. As in the bounded domain case, E is a right inverse of r^. Then we may specify 
g E H^/^(7) anduy := Eg E H^i{il) as well as our variational formulation for u = u + Eg: 

o 

Find u E V\\^(Q), such that 

B{u, w) := {AVu, Vw)^ = if, w)^ - {AVEg, Vw)^ =: F{w) ^wE Hi,(fi). 
Finally, we introduce 

D{Q) := {if E I divif E , 

which is a Hilbert space with respect to the canonical scalar product 
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2 Upper bounds for the deviation from the exact so- 
lution in dimensions > 3 

o 

Let V be an approximation of m G H^]^(f2) + u-y C Hi^(f2), where v is assumed just to 
belong to Hi^(f2) since the boundary condition may not be satisfied exactly. Our goal is 
to obtain upper bounds for the difference between Vu and Vf in terms of the norm 

MA,n ■= l^^^VIL = {A<^,'^)T ■ 



We use ( [L4| ) and get for all w e H'i^{^l) 

{AV{u - v), Vw)^ = (/, w)^ - {AVv, Vw)^ . (2.1) 
Before we proceed we note two useful results. 
Theorem 3 Let u,v G H|'i-^(i7) be as above. Moreover, let ^ be a linear and continuous 

o o 

functional over H^^^Q) and 0$ > 0, such that for all w G Hi.;^(fi) 

{AViu - v), \/w)^ = ^w) < c$ II VwlU^f^ 

holds. Then 

\\V{u - v)U^^ < c$ + 2 ||V(w - v)U^^ (2.2) 
for aim, V G Hi]^(r2), for whichu—v coincides withu—v on the boundary-/. If additionally 

o 

u — V belongs to H^^{Q) then 

\\V{u-v)U^^<c^. (2.3) 

Proof We consider 

o 

w := u — V — {u — v) G Hi]^(f2). 
Using Cauchy-Schwarz' inequality we obtain 

II Vu'i^,^ = {AViu - v), Vw)^ - {AV{u - v), Vw)^ 
< (c^ + ||V(n - v)\\A,n) WMA,n 



and thus ||Viy||^Q < c$ + ||V(u — ^')IUq- By the triangle inequality we get (2.2). (2.3) is 



trivial since we can set w := u — v, i.e. u := v := 0. □ 
We may be more specific using the trace and extension operators from Remark [2j 
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Corollary 4 Let the assumptions of Theorem be satisfied. Th 



en 



- v)\\^^^ < c$ + 2 \\VE{g - T^v)\\^^^ < c$ + 2c^ \\g - t^v\\^. 
Here > is the constant in the inequality 

||Vi5<^IU,f,<c^||</^||HV2(^) V^gH^/2(7). 



(2.4) 



Proof Setting u := Eg and v := Er^v as well as using (2.4) proves the inequalities. We 
note that (2.3) follows directly from the corollary as well. □ 



In the subsequent sections we introduce and discuss some different functionals $ and 
corresponding constants c$. 



2.1 First estimate 

o 

For any y G D{Q) and any w G Hi-^(i7) we have 

{div y,w)^ + {y,Vw)^ 



0. 



Combining (2.1) and (2.5) we obtain for all w G Hi]^(i7) and all y G D{Q) 

{AV{u-v),Vw)^ = {f + dwy,w)^ + {y - AVv,Vw)^ =: ^{w). 



(2.5) 



(2.6) 



By Cauchy-Schwarz' inequality, (1.3) with '■= 2/(A^ — 2) and (1.5) we estimate the 
right hand side ^{w) of (2.6) as follows: 

\{f + divy,w)^\ < ||/ + div?/|L ^ M-i.n < ll/ + divy||io WMln (2-7) 



< 



Cat 



\{y - AVv,Vw)J < \\y-AVv 
By Corollary |4] we arrive at the following result. 
Proposition 5 Letu,v be as in Theorem^ Then 



(2.8) 



||V(M-t;)|Uo< 



Cat 



11/ + divyli^f^ + \\y - AVw||^_i_f^ + 2c^ \\g - t^v\\^^ 



/2 



(7) 



(2.9) 



where y is an arbitrary vector field in D(f2). 

Remark 6 If v satisfies the prescribed boundary condition, then (2.9) implies 

On 



\\f + diYyl^ + \\y-AVv\\ 



n ■ 



(2.10) 



The estimates (2.9) and (2.10)) show that deviations from exact solutions of exterior 
boundary value problems have the same structure as for problems in bounded domains, 
namely they contain weighted residuals of basic relations with weights given by constants 
in the corresponding embedding inequalities. 
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2.2 Second estimate 



Assume that f2 is decomposed into two subdomains f2j and f2e with interface T := d^lg 
(see Figure [T]) and that the fields y E D (il) exactly satisfy the relation 



divy + / = in Q^- 



(2.11) 



In particular, this situation may arise if the source term / has compact support and y is 
represented (in the exterior domain Q^) as a linear combination of solenoidal fields having 
proper decay at infinity. In this case, the estimate of Proposition [5] turns trivially to 



II - ^)IU,n < Co 1/ + divyl^^ + \\y - AVvW^.i^^ + 2c^ \\g - t^v\\^,^^ 



(7) 



(2.12) 



which holds for all y G D(f2) additionally satisfying (2.11), where the weight constant is 

C7v(l + ||r|Lo,) 



Co : = 



(2.13) 



which follows directly from 

11/ + divyii^j^ = 1/ + div?/||^_^^ < IpIooa 11/ + divylj^^ < (1 + \r\oo,nj \\f + divyW^^ . 

But we also may derive another estimate. We rewrite ( |2.7 ) and use Cauchy-Schwarz' 
inequality in f2j 



\{f + divy,w)^\ = \{f + dwy,w)^^\ < \\f + divy\\^^ 



\w\ 



and estimate 



\\Vw\ 



(2.14) 



(2.15) 



Here cq. denotes a Poincare/Friedrich constant associated with the bounded domain Qi, 
i.e. the best constant of the inequality 

Mn, < cn, WVifW^^ V(/^ G 1^ G H\Qi) \ TgnMi = on 7} , 

where tqq^ : H^(fii) — )■ H^^'^(d^li) denotes the trace operator. In this case, we have again 
(2.12) but now with the (optional) weight constant 

" — (2.16) 



We note that the constant (2.13) may also be achieved by (2.7) and the argument (2.14) 
if we replace the estimate (2.15) by 



Ikin, < (1 + klooA) IklLiA < (1 + klooA) \\M\-i,n < + I'"! 



ca 



llVwl 



A,n ■ 



We summarize and get our second a posteriori error estimate. 
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Proposition 7 For ally G D{Q) with (2.11) we have 



- ^)IU,Q < Co 1/ + diYyj^^ + \\y - AVvj^^r^^ + 2c^ \\g - r. 



7^IIh1/2(-^) ) 



where Co is defined either by (2.13) or by (2.16) 



Remark 8 In general, the number cq. will be smaller and thus provides a better bound 
than cn{1 + On the other hand, the number cn{1 + Ikioo qJ/v^ ^■^ csn easily 

computable upper bound for the best possible constant Co- 



2.3 Third estimate 

Let yi and ye be the restrictions of some y G L^(f2) to fij and fig, respectively. Assuming 
yi G D(nj) and ye G D(ne) but not necessarily y G D{Q) we use the equations 



{yi,Vw)^^ + {div yi,w)^^ = {Tn,ryi,rrw)^ , 
{ye, Vw)^^ + (div ye, w)^^ = - (r„,rZ/e, rrw)^ , 



(2.17) 
(2.18) 



which hold for all w G Hi]^(^7) and in the sense of the traces rp : H^^(r2) — t- H^/^(r) and 
Tn,r ■ D(fii) H"^/^(r) respectively r„,r : D(fie) H"^/^(r). At this point we assume 
that the interface F is Lipschitz (in order to guarantee that the traces are well defined). 
By {Lp,-ip)^ we denote the duality product of H"^/^(r) and H^/^(r). We recall that the 
normal traces r„,r2/i and Tn^rVe possess weak surface divergences in H~^/^(r) as well. If 
y G D(n), then divyi = divy in fij and div ye = divy in fie- Hence, in this case adding 



(2.17) and (2.18) we obtain by (2.5) 



{Tn,ryi - Tn,rye, Trw)^ = {y, Vw)^^ + {divy, w)^ = 

o 

for all w G Hi]^(f2). Therefore, we get 

rn,ryi = Tn,rye 

for all y G D{Q) since rr is surjective. 



On our way to find $ like in (2.6) we now insert (2.17), (2.18) instead of (2.5) into 



(2.1 ) and obtain 



{AV{u-v),Vw)^ = {f + diYyi,w)^^ + {f + diYye,w)^^ (2.19) 

+ {y- AVv, Vw)^ + {TnxVe - rn,ryi, Trw)^ =: $(w). 



The third term of $(w) will be estimated by (2.8) and for the last term we may use 
the continuity of the trace operator rr in combination with a Poincare/Friedrich estimate, 
i.e. 



o 



(2.20) 
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and obtain 



{'Tn,rye ~ Tn,r2/i) ''"7w)p| < iTn^rVe ^ Tn,rZ/i || h-i/2 (r) II ''"r'?^ || hi/2 (r) 

< cr ||r„,r?/e - r„T2/i|lH-i/2(r) l|Vw" 



(2.21) 



A,n ' 



To estimate the second term of ^{w) we again use (1.3) and (1.5) and obtain 



|(/ + div|/e,w)o I < ||/ + div?/e|li. 



o IklLio < ||/ + div?/e||io |k||_io 



(2.22) 



< 



Cn 



||/ + divi/e|L o ||Vw| 



An 



Considering the first (and last) term of $(w) we have once more at least two options 
as in section 12.21 to obtain the estimate 



|(/ + div?/i,u;)^J < Co \\f + divy,\\^^ WM\A,n 



(2.23) 



with Co defined either by (2.13) or (2.16). 

Finally with ( [2l9| and (|2.8D, ( |2!2Tf , ( |z!22| , ( [2^ we get by Corollary |4| the third 
estimate. 



Proposition 9 For all y G L^(i7) with yi G D(f2j) anc? ye G D(f2e 

Cat 



l|V(^^-t^)IUo < Co||/ + div?/i||o. + — = ||/ + div?/e||i 



\\y 



we have 



(2.24) 



+ Cr ||Tn,rl/e ~ Tn,r2/i|lH-i/2(r) + ||(7 — T^f |||_|i/2^^-) 

wt/i Co /rom Proposition T/ie ng'/ii /ianti szc^e of (2.24) vanishes if and only if v 
coincides with u and y with AWu. 



Remark 10 There are many ways to deduce (2.20). We just mention that Tr<f can be 
considered as a trace of a function defined in Qi or or even of a function, which is just 
defined in a small neighborhood ofT. Thus, we may adjust the constant cr according to 
our needs. 

Remark 11 This estimate suggests even a solution method: We construct approxima- 
tions using locally supported trial functions in Qi, e.g. FEM, and utilize global approxima- 
tions properly behaving at infinity for Q^- These two types of approximations are usually 
difficult to meet together exactly on the artificial boundary T. However, Proposition 
shows that this is not required because we can use instead the penalty term with known 
penalty factor cr- In addition, we have one more parameter, the 'radius' of the interface 
r. Since T is artificial and arbitrary we can use this parameter in the algorithm in order 
to obtain better results. 

Remark 12 At this point we shall note that all our estimates are sharp, which easily 
can be seen by setting v := u E and y := AVu G □((]). 

Remark 13 In Propositions\^ \^\^we can always replace the last summand of the right 
hand side by 2 ||V('U — "0)1^^ or 2 \\'VE{g — Tyv)\\^^ using Theorem^ and Corollary^ 
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3 Upper bounds in dimension N = 2 

Of course, Theorem [3] holds for = 2 as well and the modifications on the estimates 
depend just on the Poincare/Friedrich estimate and thus they are obvious using the proper 
Cauchy-Schwarz inequality. We achieve 

Proposition 14 Let C M^, such that M.^ \Q contains the unit ball. 
(i) For all y G D{n) 

2 

- ^)IU,n < \\rlnr{f + dwy)\\^ + \\y - AVv\\^^,^^ + 2c^ ||^ - t^vW^,^^^^^ . 



(ii) For all y G D{Q) with divy + f = in Qf. 

- ^^)IU,f7 < Co 11/ + divyllf^^ + \\y - AVv\\^-,^^ + 2c^ \\g - t^v\\^,^^^^^ 

where Co = min |2 ||r lnr||^ , cn,| /^/ca- 



(iii) For all y G L (fi) with yi G D(f2j) and ye G D(f2e) 

2 

l|V(M-f)IU,n < Co \\ f + div yi\\^^ + -— ||rlnr(/ + divye)||i,n, + lb - ^Vt;|U_i 



+ Cr \\Tn,rye — Tra,rl/j|lH-i/2(r) + \\g — r^f |||_|i/2(^-) • 



Analogously, Remarks]^ 10, 11, 12, 13 hold. 



A Appendix 

A.l Lower bounds for the error 

We note by a standard variational argument 



l|V(w-t^)|lln= sup (2{AW{u-v),y), 



Thus, we obtain for all w G H|':]^(i7) the estimate 

||V(n - v)\\X^ > 2 {AV{u - v), Vw)^ - \\Vw\\X^ 

= 2 {AVu, Vw)^ - {AV{2v + w), Vw)^ , 

which is sharp since one can put w = u — v. But to exclude the unknown exact solution 



u from the right hand side we need w G H^_^(fi) since then by (1.4) 



I V(« - v) 111 o > 2 (/, w)^ - {A\/i2v + w), Vw)^ . (A.l) 
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But this estimate is no longer sharp because we can not put w = u — v anymore. In fact, 
with AVu e D(fi) and div AVu = -f we get for w G H^i(i7) 

{AVu, Vw)^ = (/, w)^ + {Tn,^AVu, T^w)^ . 

Hence, we obtain the estimate 

II V(m - v) > 2 (/, w)^ - {AVi2v + w), Vw)^ + 2 (r„,,AVM, r^w)^ 



for all w G H^^(f2), which is sharp and coincides with (A.l) if w G H^]^(fi). But the 
unknown exact solution u still appears on the right-hand side, i.e. the normal trace of 



AWu on 7. Furthermore, if (r„^^y4Vu, r^w) > then (A.l) can not be sharp 



A. 2 Poincare type estimates for exterior domains 

We introduce the radial derivative dr := ^- V, where ^(x) := x/r{x). Furthermore, -B^ and 
Se denote the open ball and sphere of radius e centered at the origin in M^, respectively. 
We will use the ideas of [H Lemma 4.1] and fH Poincare's estimate III, p. 57] with some 
minor useful modifications. 



Lemma 15 Let Vt d , N > I, he a domain and /3 G M. For all u G C°°(fi) the 
following Poincare estimates hold: 



(i) If 13 > I- N/2 then 



{213 + N - 2) ||r'^"^M|L < 2 llr^ dr 



u\ 



(ii) Let Bi(im.^\VL. If (3 > {?, - N)/2 or (3 <l - N/2 then 

.13-1 



\2l3 + N-3\ 



Inr 



u 



< 2 \\r^ drU 



(iii) If N = 1 then 



|2/3- 1| ||(l + r)^"^M|| < 2||(1 + r)''9^M|| + |2min{0,2/3- 1}|'^^|m 



1/2 



where u will he extended hy zero to M. 



For the estimates derived in this paper it suffices to set /3 = 0. In this particular case, 
the above lemma implies 

o 

Corollary 16 Let Vl C M^, N > 1, he a domain. For all u G C°°(fi) the following 
Poincare estimates hold: 
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(i) IfN>3 then 

2 2 

(ii) If N ^2 and Bi cM.'^\n then 

\\u/{rlnr)\\^<2\\dru\\^<2\\Vu\\^. 

(iii) IfN^l then 

||«||_i,a < ll«/(l + r)\\a < 2 \\dru\\^ + ^/2\u{0)\ < 2 \\u'\\ 
Hence, if C R± we have 

\\uU^^<\\u/{l + r)\\^<2\\dru\\^<2\\uX. 

Remcirk 17 Of course, by continuity all these estimates extend to appropriate weighted 
-Soholev spaces. 

o 

Proof Let Q. C M^, iV > 1, be a domain and u e C°°(n). By partial integration we get 
for all a e M and £ > 



2 / r'^udrUdX^ I r'^ dr\u\^ 
Jn\Bs Jn\Be 

-{a + N-1) [ r"->|MA-£"/" \u\'^ da. 

Jn\Be JSe 



Thus, for all 7 e M and j3 :^ {a + l)/2 



r^'u dr u dA 



u\ 



|2 

il\B, 



+ 7(7 - 2/3 - iV + 2) ||r^-'«||^\5^ - ie^^~^ J \u\' da. 



Now the left hand side of this equality converges by the monotone convergence theorem. 
Since r*^ G L^(f/i), if and only if z/ > — A^, and | {u]"^ da\ < ce^~^ the right hand side 
converges for /3 > 1 — N/2 by Lebesgue's dominated convergence theorem in M. Hence, 
for £ ^ we obtain 

\y dru + -fr^~^u\\l = \y dru\\l + 7(7 - 2/3 - + 2) ||r^~'ii|£ ■ 
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Choosing ^ :— 2/3 + N — 2 > we finally get by the triangle inequality 

Since we are especially interested in the case (3 = this estimate is only applicable in 
dimensions N >3. 

For = 1 we proceed as follows: For all a e R we have 

2 [ {l + r)'^udrudX = ±2 I {1 ±tYu{t)u{t)' dt 

= ±2 I {l±tY{\u{t)\^)' dt = -a I {l±tY-^\u{t)\^dt-\u{0)\^ 
Jr± Jr± 

and thus 

2 [ {l + r^udrUdX^ -a [ {1 + ry-^\u{t)\'^ dX - 2\u{0)\'^. 
Jr Jr. 

Hence, for all 7 e K and /3 := (a + l)/2 



|2 

In 



1 + rf dr u\\^ + |7|2 II (1 + r)f-\\\^ + 27 ((1 + rf dr u, (1 + rf-'u)^ 

" V ' 

= [ {l + r)''udrudX 
Jn 

= 1(1 + rf dr u\\l + 7(7 -2^ + 1) 1(1 + rf-'u\\l - 2j\um'. 

As before the triangle inequality and the choice 7 := 2/3 — 1, but now without any 
restrictions on (3, lead to 

I7I l(l + r)'^-^M||^ < ||(l + r)^9,M||^+ (\\{l + rydru\\l-2^f\u{0)\''y^^ 
< 2l(l + r)'^a^'ul^ + |2min{0,7}|^/^|'u(0)|. 

The remaining case N = 2 requires the use of logarithms. Moreover, the origin is now 
a problematic singularity, which has to be removed from our domain. Therefore, we may 
assume Bi C R-^\Q and N > 1, having A?" = 2 in mind. We start once more for all a e R 
with 

2 / - — udrudX= / - — dr\u\^dX 

= -(a + N-l) / \u\^dX + / ^kPdA. 

' Inr ' ' Ja InV' 
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Now our usual procedure gives for 7 e R and ^ := (a + l)/2 > 



drU + 7- u 

mr 



= ||r^9r it||f^ + I7 



Inr 



u 



+ 27(r^ drU, 



Inr 



u 



= llr^^rMlL +7(7 + 1) 



,,3-1 



Inr 



-u 



- — u dr u dA 
mr 



-i{N + 2P-2) 



.13-1 

/inr 



zU 



Thus, for 7(A'" + 2^ — 2) > we can estimate 



drU + ^ 



Inr 



u 



< ||r^(9,M||' +7(7-2/3- AT + 3) 



Inr 



u 



which leads to the estimate 



Inr 



if we set 7 := 2/3 + — 3 with the additional constraint 7(7 + 1) > 0, i.e. 7 > or 
7 < — 1. Finally, again by the triangle inequality 



|7| 



.13-1 



Inr 



u 



< 2 ||r''9^M||fj 



follows for aU /3 > (3 - A^)/2 or /3 < (2 - N)/2. 



□ 
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